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GRAPH HOMOLOGY AND STABILITY OF COUPLED OSCILLATOR NETWORKS 

JARED C. BRONSKI , LEE DEVILLE , AND TIMOTHY FERGUSON* 

Abstract. There are a number of models of coupled oscillator networks where the question of the stability of fixed points 
reduces to calculating the index of a graph Laplacian. Some examples of such models include the Kuramoto and Kuramoto— 
Sakaguchi equations as well as the swing equations, which govern the behavior of generators coupled in an electrical network. 
We show that the index calculation can be related to a dual calculation which is done on the first homology group of the graph, 
rather than the vertex space. We also show that this representation is computationally attractive for relatively sparse graphs, 
where the dimension of the first homology group is low, as is true in many applications. We also give explicit formulae for the 
dimension of the unstable manifold to a phase-locked solution for graphs containing one or two loops. As an application, we 
present some novel results for the Kuramoto model defined on a ring and compute the longest possible edge length for a stable 
solution. 


1. Introduction. Let G = (V^E^T = {7i,,,u}i,,ujgu) be a weighted graph. The Laplacian matri3|^of G 
is the \V\ X \V\ matrix whose components are 


i^o) 


VW 


V ^W, 
-Hulvu, V = w. 


A classical result is that if 7^,^, > 0, then Cq is negative semidefinite, with the number of zero eigenvalues 
equal to the number of connected components of the graph G. If some of the 7„u) are negative, then Cq can 
have positive eigenvalues as well. In [^[^, the authors consider the problem where some subset of the 'yvw 
are negative, and show upper and lower bounds on the number of positive eigenvalues of Cq that depend 
only on the “sign topology” of the graph, i.e. the pattern of which edges are positive and which are negative, 
regardless of their weights. 

More specifically, the main result of is: assume that the graph is symmetric (7(„,m) = 7(ti,,u)), and 
define G+ (resp. G_) to be the subgraphs of G containing only edges with positive (resp. negative) weights. 
Let |G| be the number of components of a graph G. If n+(G) is the number of positive eigenvalues of Cq, 
then 


|G+|-l<n+(G)< |U|-|G_|. 


( 1 . 1 ) 


Moreover, one can show that these bounds are saturated, i.e. for any integer in that range, there is some 
choice of weights that give exactly that many positive eigenvalues. (In fact, one can say much more about 
the genericity of the sets of weights that give these values, see for more detail.) 

It follows from the formula above that if G is a tree, then the number of positive eigenvalues is exactly 
the number of negatively weighted edges in the graph. To see this, notice that, for a tree, the number of 
connected components of G+ is exactly one plus the number of negative edges: if we start with a tree with 
only positive edges, every flip of an edge from positive to negative will add another component to the positive 
subtree. Conversely, since the negative edges cannot form a cycle, the number of connected components of 
G_ is the number of vertices minus the number of negatively weighted edges. In short, if G is a tree, then 
n+(G) is independent of the weights of the edges, and very easy to compute by inspection. 

Now, consider the case where the graph has one cycle and one negative edge. It is easy to see in this 


case that the upper bound in (1.1) is one, and the lower bound is either zero (if the negative edge lies in 


the cycle); or one (if the negative edge does not lie in the cycle). In the latter case, the number of positive 
eigenvalues is again independent of the weights. In the former case, the weight on the edge matters, and it 
is not hard to imagine that the relative magnitude of the weight on the negative edge compared to the other 
edges in the cycle determines the stability. 

From these examples, it is at least plausible to conjecture that the stability computation is something 
that can done on a cycle-by-cycle basis. This is further supported by previous work of the authors [^, where 
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it was shown that a certain quotient of Hi(G) by a subgroup was important for understanding the spectrum 
of the graph Laplacian with positive and negative edges. 

In fact, this is the result of the current paper. Let us denote the number of cycles by C = \E\ — \V\ + 1. 
Then the main result of this paper is that the computation of the index of Cq is the equivalent to computing 
the index of a particular C x C matrix (what we call the “cycle intersection matrix”). If C <C \V\, then this 
is clearly a much simpler problem; the technique we describe here is most effective for sparse graphs. 

1.1. Connection to nonlinearly coupled phase oscillators. As stated above, the problem we 
consider is a pure linear algebra problem and can be stated as such, but in fact it is inspired and informed 
by considering its “killer app”, i.e. understanding the dimension of the unstable manifold of a fixed point of 
a nonlinearly coupled network dynamical system. 

Some examples of the kinds of models we are interested in include the Kuramoto problem on a general 
weighted connected graph D = {Sij}, namely: 


d0. 


^ 6ij sm{9j -9i). 


( 1 . 2 ) 


dt 

More generally one can also consider the Kuramoto-Sakaguchi model 

=u},+^5ijsuv{9j-9i + a), (1.3) 

as well as the swing equations, which govern the dynamics of a coupled system of generators and loads 


+ XI sin(% - 9i). 


(1.4) 


In each of these examples the dimension of the unstable manifold to a fixed point of the equations is given 
by the number of positive eigenvalues of a graph Laplacian. 

The study of the fixed points of these systems, and their stability, has lead to a huge number of papers 
in the literature [6 -23 . It has been widely observed that the behavior of solutions varies dramatically as 
the topology of the graph changes: these models are essentially trivial on a tree, but become more complex 
as the graph becomes denser. This again suggests that the first homology group Hi{D), which counts the 
number of loops in the graph, should be important 

To date, most of the analytical work has focused on special graph configurations with a high degree of 
symmetry. The most-studied case, where the behavior of the model is very well understood, is the case of 
“all-to-all” coupling, where the underlying graph is the complete graph , and all of the edge weights Stj are 
equal. In a ground breaking paper 19 , Mirollo and Strogatz gave upper and lower bounds (which differed 


by 1) on the dimension of the unstable manifold to an arbitrary fixed point solution to this model. Dorfler 
and Bullo 24 gave a sufficient condition for the existence of a stable phase-locked solution, and two of the 


current authors 25 gave an explicit formula for the dimension of the unstable manifold. Similarly Verwoerd 
and Mason give a proof of the existence of a critical coupling constant for the cases of the complete 16 and 
bipartite graphs, and one of the current authors [26] gives some results on the minimizers and saddles for 


the Kuramoto on cyclic graphs. For general topologies, it was first observed by Ermentrout 27 that if the 
graph is connected and all of the links are short (the angle differences are less than |) then the fixed point 
is stable. 

We show in this paper that the intuition put forward in the preceeding paragraphs is correct. We show 
that, rather than computing the number of positive eigenvalues of the Laplacian, an object that acts on 
the vertex space of the graph, there is a dual calculation that can be posed on the cycle space. This dual 
problem gives an alternative method for computing the dimension of the unstable manifold to a particular 
fixed point. This approach is particularly attractive for problems in which the graph is relatively sparse, 
and the dimension of the cycle space is much smaller than the dimension of the vertex space. We note that 
there are some very important applications for which this condition is satisfied. One example is the graphs 


associated with electrical power networks. Such networks are governed by the swing equations (1.4) and these 
tend to be reasonably sparse. The following data is from Cotilla-Sanchez, 


Hines, Barrows and Blumsack , 
who calculated the number of edges and vertices for the IEEE 300 model (a standard test power system) as 
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well as the Eastern (El), Western (WI) and Texas (TI) interconnects, the largest components of the North 
American electrical network^ 


Table 1.1 

Sparsities of Various Electrical Networks 


Network 

Edges E 

Vertices N 

Number of Cycles C = E — N + 1 

Sparsity % 

IEEE 300 

409 

300 

no 

.37 

El 

52075 

41228 

10848 

.26 

WI 

13734 

11432 

2303 

.2 

TI 

5532 

4513 

1020 

.23 


As can be seen from the data above electrical networks tend to be quite sparse, with the dimension of the 
cycle space being 20—40% of dimension of the vertex space. Eor instance if one wanted to understand the 
stability of a fixed point of the Texas Interconnect, the dual problem presented in this paper would reduce 
the calculation to the eigenvalues of a 1020 x 1020 matrix, rather than the 4513 x 4513 eigenvalue problem 
determined by a naive stability calculation. 

2. Summary of Results. In this section, we define our notation and state the two main theorems of 
the paper, and give a sense of the main ideas of the proofs. The complete proof of the main results will be 
deferred to Section [5l 

2.1. Definitions. Definition 2.1. A weighted graph G is the triple G = {V,E,T), where V is a 
set of vertices, E C V x V a set of edges, and T = {yejegB a set of edge weights. By this definition, the 
graphs in this paper are always simple (at most one edge between vertices) and undirected. 

Definition 2.2. Let G = {V,E,T) be a weighted graph. The Laplacian matrix of G is the \V\ x \V\ 
matrix whose entries are given by 




V,W 


1 


V w, 

V = w. 


Since the graph is undirected, the Laplacian matrix is always symmetric and has real spectrum. We denote by 
n+(G), no(G), n_(G) the number of positive, zero, and negative eigenvalues of Lq, and refer to this triple as 
the spectral index ofG. Since, for any G, CqI = 0, Cg has a zero eigenvalue and thus a zero determinant. 
We define the reduced determinant of Cq, det|.ed(GG), to be the determinant of the matrix restricted to 
the space of zero sum vectors, i.e. to the orthogonal complement of \ in (Note that det|.ed(>C( 3 ) = 0 iff 

Cg has a nonsimple eigenvalue at 0.^ We will colloquially refer to a Laplacian as stable if n+{CG) = 0 and 
unstable otherwise. 

Definition 2.3. Let G = G{V, E,T) be a connected weighted graph. Choose and fix and orientation for 
each edge (so that each edge has both a head and tail vertex). Define the signed incidence matrix of G 
as the \V\ x \E\ matrix whose entries are given by 


{ 1, e = (u,*), 
- 1 , e = {*,v), 
0, else. 


The nullspace of the matrix is called the cycle space. It is a standard result l2tfj that if G is connected, 
then the dimension of this nullspace is C := |£’(G)| — \ V (G)| + 1. We call this integer the cycle rank of the 
graph. Note also that while the definition of B depends on the choice of orientation, the number C does not. 

Using the Rank-Nullity Theorem, the above is equivalent to saying that B^ has nullity one, a fact we also 
use below. It is also not hard to see that B^l = 0, so in fact kei{BQ) = span(l) whenever G is connected. 


small amount of processing has been done on this data, including combining parallel edges and removing disconnected 
buses. Other sources give slightly different numbers without changing the general picture. 
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We define the tree set Tree(G) to be the set of edges of G that are contained in every spanning tree 
of G. The complement o/Tree(G) is the cycle set Cycle(G), consisting of every edge which is included in 
some cycle in the graph. In terms of the incidence matrix, e is in the cycle set if and only if there is a vector 
X S N{Bq) with Xg 7^ 0. 

Definition 2.4. Let yi, 1 / 2 , ■■ ■ ,yc be a basis for the cycle space. Define the \E\ x G matrix Yc whose 
columns are given by the yi, and Dq to be the \E\ x \E\ diagonal matrix where the eth entry is 7 e. Finally, 
we define the cycle intersection matrix, or the cycle form 

ZG-.= -Y(fDa^YG, ( 2 . 1 ) 


i.e. the matrix with components 


(Zah = - XI 

eeE 


yi,eyj,e- 


Note that Zg is always a C x C matrix. The exact form of Zg will depend on the choice of basis, and we 


will specify the exact choice of basis in (5.3) below. For now, we simply note that any property that we need 
for Zg will be independent of this choice. 


Lemma 2.5. Let G = (V,E,T) be a connected weighted graph. Define D as in Definition 2.4 and B as 
in Definition\2.3\ Then 


Lg = —BgDgBq. 


Proof. Let us first assume that v ^ w. Then 

{BgDgBg)v,w = X BygDee'Bji^g = ^^BveBwele- 

e,e' e 

If v,w are adjacent, then either e = {v,w) or e = {w,v) is an oriented edge. In either case, By^B^^g = —1. 
Moreover, if e is an edge not incident on both v and w, then BygBy,g = 0, and thus the sum has a single 
term, and {BGDGB]^)yy, = —^g. Therefore all of the off-diagonal terms are the same as —Cg- 

Finally, notice that since B]jl = 0, we have BgDgBqI = 0, and from this it follows that the diagonal 
terms must be correct as well. 

□ 

2.2. Theorems. The following are the main results of the paper. The first result gives a formula for 
the number of positive eigenvalues of the Laplacian in terms of the number of edges with negative edge 
weights together with the number of positive eigenvalues of the cycle intersection matrix. 

Theorem 2.6. Let G = {V,E,T) be a connected, weighted graph and let Cg be its Laplacian matrix. 
Then the number of positive eigenvalues of the Laplacian is the number of negative edges in the graph minus 
the number of positive eigenvalues of the cycle intersection matrix, i.e. 

n+{CG) = MeCE\-fg<G}-n+{ZG). ( 2 . 2 ) 

Recall in the definition of Zg that one makes a choice of basis; however, it follows from Lemma \A .3\ that the 
index of Zg is independent of this choice of basis. 

In particular, we have the inequality 


n+{CG) < #{e e Tree(G) | 7 e < 0}. 


From (2.2), the maximal number of negative edges a graph can have and still have a stable Laplacian is the 
cycle number G. In the case where the Laplacian is stable and there are exactly C negative edges, then there 
is choice of basis cycles so that there is exactly one negative edge in each cycle in the basis. In other words 
if there are exactly G negative edges, then the graph with the negative edges removed must be a tree. 

Remark 2.7. The basis for the cycle space can be chosen in the following way: Let G+ be the subgraph 
containing only the edges of G with positive weight. The cycle space of G+ forms a subspace S'+ of the full 
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cycle space of G. Similarly let G- be the subgraph containing only negatively weighted edges, and S- the 
corresponding cycle space. Note that and S- are necessarily orthogonal. The basis for the cycle space of 
the full graph G can be chosen as the union of bases for 5'+ and S- together with a set of “mixed” cycles 
that cannot be decomposed into a sum of cycles containing only edges of one sign. These mixed cycles are 
a basis for the quotient of the full group of cycles by the subgroups S±. In this basis the cycle intersection 
form can be written in the following partitioned form 


Am 

Bm,+ 

Bm.- 

^M,+ 


0 

Bit,- 

0 



Zg = 


where and Zq_ are the cycle intersection forms for the subgraphs G.^. and G- respectively. 

First we note that Zq^ is negative definite and Zg_ is positive definite, so that n+^Zo) is at least as 
large as the dimension of the cycle .space of G- and no larger than the dimension of the full cycle space 
minus the dimension of the cycle space of G+. These inequalities are equivalent to the inequalities derived 
previously by the authors 0 

Secondly we note that the calculation of the index of Zq can be further reduced to a calculation on the 


space of mixed cycles. In particular by applying the Haynsworth theorem to (2.7) we have that 




and thus the calculation can be reduced to finding the index of a matrix whose .size is the number of mixed 
cycles. However since the above formula requires inversion of the matrices Zq^ it is not obvious that this is 
more attractive than working on the full cycle space. 

The second theorem is a similar result at the level of determinants rather than the spectral index. This 
result is not really needed for the applications, but it is proven with similar ideas and provides a nice dual 
form of the famous matrix tree theorem. This is essentially a weighted version of a theorem proven by 
Sjogren, 


30 though the proof given here is quite different. 


Theorem 2.8. Let G = (y,i?,r) be connected, and define Cq as above. Let Ng = |y(G)|, and then 


W ^ det(Z(G)) X 7 e 

^ e^E 


(2.3) 


The quantity det(2^(G)) is a homogeneous Laurent polynomial in the edge weights which generates labelled 
tree complements: each term in the Laurent polynomial represents a set of edges that, when removed from 
the graph, yield a tree. 

Remark 2.9. There is an obvious connection between the two theorems, but they are not quite equivalent 
since they address different quantities. For example. Theorem \ 2.^ seems stronger since it gives an actual 
number, but it only implies Theorem \2.6\ modulo 2. 

We defer the proof of these results until a later section, but we will briefly summarize the two main ideas 
here. The first is that for any G, there is a covering tree T (which is a weighted tree), and a graph map 
Lp: T —>■ G, which is surjective on vertices and bijective on edges. In other words if one identifies certain 
vertices of the tree T one obtains G. The identification of vertices in the graph corresponds at the level of 
linear algebra to restriction of the Laplacian to a certain subspace, so the Laplacian on any graph can be 
realized as the restriction of the Laplacian on a tree to an appropriate subspace. 


The second idea is to use the Haynsworth theorem 31 , which is stated in Lemma A.l The Haynsworth 


theorem is a kind of duality result relating the number of positive eigenvalues of a matrix to the number 
of positive eigenvalues oi a, k x k principle minor matrix and the number of positive eigenvalues of the 
complementary {n — k) x {n — k) minor of the inverse matrix. In other words, modulo some invertibility 
assumptions, if M is self-adjoint and S' is a subspace of then 

n+{M) = n+{M\s) +n+(M~^\g±). 

Brushing aside some technical points associated with invertibility which will be addressed below, the basic 
observation is that if we take M to be the Laplacian on the covering tree and M\s to be the Laplacian on the 
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graph G, then one can compute that M~^\g± is exactly the cycle intersection form Zq- For a tree it is easy 
to see that the number of positive eigenvalues is exactly the number of edges with negative edge-weights, 
and formula (2.2) follows. There is an analogous identity at the level of determinants, giving (2.3). 

3. Linear Examples. 

3.1. Diamond graph. We consider the simplest graph with cycle rank (7 = 2, the diamond graph 
shown in the left frame of Figure [5T| with five edges and four vertices. The Laplace matrix is given by 


C-G = 


Id 

0 




( -{la+ld+le) la le 

la -{la+lb) lb 

le lb -{ib + lc+le) 1, 

\ Id 0 7c -{ic + ld) J 


By the matrix-tree theorem, we can obtain the eight labelled spanning trees by considering the 4,4 cofactor 
det£44((7) = 


-(7a+7d+7e) la le 

la -ila+lb) lb 

le lb -{ib+lc + le) 

= lalblc + lalbld + lalcld + lalcle + laldle + Iblcld + Iblde + Ibldle- 


(3.1) 


Let us orient this graph clockwise around the outside, then 1 - 
cycle space is given by 

zi = (1,1,0,0,-1)‘ 

Z2 = (0,0,l,l,l)‘. 


3. With this orientation, a basis for the 


The cycle intersection matrix Zq is given by 


/ 


^G = - 


1 

la 


1 

lb 

1 

le 


1 

le 


1 

7c 


1 

le 

1 1 

-1- 

Id le 


The determinant of this matrix is 
det Zq = ^ 

Idle 


Icle 


Ible 


Ibid Ible lale laid laic 


(3.2) 


It is not hard to see that if we multiply (3.2) by lalblcldle, then we obtain the expression in (3.1), verifying 
Theorem 2.8 The Aei{ZQ) generates labelled tree complements: each term in this Laurent polynomial 
represents a pair of edges that, when removed from the graph, yield a tree. 

Moreover, we can use Theorem |2.6| to compute bounds on the edge weights to determine stability. For 
ease of notation, we write pi = l/ 7 i, and 


Wl=Pa+Pb, W2=Pc + Pd, Wi2=Pe- 

Then we have 

^ ^ _ f Wi+Wi2 Wi2 

Wi2 W2 + Wi2 

Let us first consider the case where there is one negative edge. In this case, Zq has either one or zero 
positive eigenvalues, which is the same as saying that —Zq has one or zero negative eigenvalues. If there is 
one, then Cq is stable, and if zero, it is unstable. 
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Note however that det Zq < 0 if and only if there is one positive and one negative eigenvalue, which 
gives the condition 


(iCl + W12)(W2 + W 12 ) — rci2 < 0, 


or 


W1W2 + W1W12 + W2W12 < 0 . 


If we assume that 7 e < 0, then, solving for W 12 , we obtain the inequality 

W 1 W 2 

W 12 < -^-, 

Wi + W2 


or 


^ Pa + Pb + Pc + Pd 
[Pa + Ph) (Pc + Pd) 

Notice, for example, that if we assume that 70 ~ 0, then the pa terms dominate, and the condition we obtain 
for stability is 7 e > —{pc + Pd)- Here 7 ^ is not relevant, and this makes sense: as 70 —>■ 0, the a, 6 path is 
almost broken as a loop and thus cannot significantly affect stability, so the condition for stability will be a 
balance between e and c, d. 

On the other hand, if one of the edges on the outside is negative (say wlog 7 a), we obtain the inequality 


Wi < 


W2W12 
W 2 W 12 ’ 


or 


{Pc Pd) Pe 
Pc + Pd + Pe 


- Pb- 


Note that every p on the right-hand side is positive. Note the particular role played by the pb term: if 
we make 7 b larger, this makes pb smaller (in absolute value), making it more difficult for pa to cause the 
instability. 

qqqqq 

Now, assume the graph has two negative edges and let us try and find conditions for stability. Note 
that we cannot have both a, b or both c, d negative, since this automatically gives an unstable Laplacian. 
Therefore we need that at least two of wi,W 2 ,wi 2 be negative. 

The condition that the matrix be stable is that —Zq have two negative eigenvalues, which we can obtain 
by the condition of having a positive determinant and a negative trace, so 


W1W2 + W1W12 W2VJ12 > 0 , Wi + W2 + ‘ 2 wi 2 < 0 . 


By symmetry, we only need to consider the case where wi,W 2 < 0 and where wi,wi 2 < 0. In the former 
case, there is a negative edge on each of the “outsides” of the graph, and in the latter there is one in the 
middle. 

3.2. One-cycle graphs. Let us now consider the general one-cycle graph and consider the condition 
for stability. Note that such a graph can have at most one negative edge and still be stable, and, moreover, 
it cannot be contained in the tree part of the graph, but must be part of the cycle. Therefore we can simplify 
and consider Rn, the ring graph on N vertices: this graph has N vertices and N edges, and let us label 
them so that edge i goes from i —b i + I, modulo N. 

The cycle can be represented by the vector 1 , and thus we have Zq is the I x 1 matrix (i.e. the scalar) 


Zg 
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Thus to obtain a stable Laplacian, we need that Zq > 0, or 


N 


.-1 


Assuming that 71 < 0 and the remainder are positive, the condition for stability is 


7i > 


-1 


n 


E 


1 ' 

7 * 


If all of the positive edges have an equal weight, 7 i = 7 , then we obtain 

7i > -l/in-l). 

In this case, we see that increasing 7 helps with stability, whereas increasing the number of edges, all things 
being equal, hurts stability. 

3.3. Two-cycle graphs. We can consider the general two-cycle graph; in fact, we can reuse many of 
the computations of Section [3d] above. 

As in the previous section, we simplify by assuming that there are no edges in Tree(G'), i.e. each edge 
is in one or more cycles. Let us split the edges into three sets S'!, 5 * 2 , 5 'i 2 ; the edges in Si are those only in 
cycle one, those in S 2 are those only in cycle two, and finally S 12 are those edges both in cycles one and two. 
It is clear that we can choose Si, S 2 to be non-empty, but for some graphs, it may then be the case that 5 'i 2 
is empty (e.g. a graph with a figure-eight structure). Let us reorder the edges so that when ordered, the sets 
Si,S 2 , S 12 come in order, and we write A:* = |S'*|. From this we see that we can choose the cycle basis as 

2/1 = ( 1 , 1 ,..., 1 , 0 , ..., 0 , 1 ,..., 1 ), 2/2 = ( 0 , 0 ,. . ., 0 , 1 ,..., 1 , 1 ,..., 1 ), 


where the three groups of edges are of length ki,k 2 , ki 2 . If we then define 


ees. 


1 

7e 


then we see that the cycle form Zq has the same structure as in (3.31, and the subsequent analysis holds 
over for this case. 

The simplest case is the one where 512 = 0 and ^12 = 0. In this case Zq is diagonal, and the number 
of positive eigenvalues is exactly the number of negative diagonal entries. For example, let us say that the 
edge e for which 7 e < 0 is in cycle 1, and thus we need u>i < 0 to have a stable Laplacian. This gives 


I 

7e 


+ E 

f^ejeSi 


1 

7/ 


< 0 , 


or 


-1 

1 /T'/ 

This is exactly the condition as in the one-cycle case. Similarly, if the graph has two negative edges, we need 
both wi,W 2 < 0 and the conditions on each cycle are decoupled. 

In the case where W 12 > 0, this gives an analysis similar to that given above. Let us consider the case 
where there is one negative edge, and it is in both cycles one and two. We again have the condition 

W1W2 

W12 < -;-. 

Wi + W2 






If, for example, we consider the case where all of the positive edges have the same weight 7 , and the one 
negative edge has weight 7 e < 0 , then we have 


wi = ki/j, W2 = ^ 2 / 7 , W 12 


ki2 - 1 J_ 
7 7e’ 


This simplifies to 


1 1 / fci 

— <-( , , , + (^12 - 1 ) 

7e 1 \ki+k 2 


(3.4) 


Again, increasing 7 helps for stability since it makes the condition easier to achieve, but increasing any of 
ki, ^ 2 , ki 2 makes the condition harder to achieve. 

The interpretation of these facts is as follows: clearly, increasing the strength of the positive couplings 


will make it harder for the negative coupling to destabilize the entire system, and in fact from (3.41 we see 


that the two weights are effectively linearly related. Moreover, we see that increasing the number of links in 
either of the cycles makes it easier to destabilize the system, similarly to the case of a single cycle. 

4. Nonlinear example — Kuramoto model on a ring. The examples in the previous section 
studied the stability of linear Laplacians, but we can also use the current results to obtain novel results for 
nonlinearly coupled oscillators as well — in particular we consider the Kuramoto model. 

It has been known for a long time that if we have a hxed point for the Kuramoto model where all of 
the links are “short”, i.e. the angle difference across a link is less than tt/2 in absolute value, then this gives 
a stable fixed point, see [^. We could then pose the question of whether there exist solutions with “long” 
links that are still stable (this would correspond to having negative off-diagonal terms in the Jacobian). This 
question was studied numerically in 21 ; there the authors considered a large class of random networks and 


random frequencies, and showed that stable fixed points with long links do exist, but they are exceedingly 
rare (typically the proportion of stable points found that have long links were less than I/IOOO). 

Intuitively, it makes sense that if we make the weight Stj on an edge arbitrarily small, then this is 
possible, but we could then pose the question: can we have a stable long-link solution if the edge weights 
are all equal? We will show that this answer is affirmative on the ring by a judicious choice of frequencies. 
In particular, we will answer the following question for the ring topology: given a specihc graph topology 
for the Kuramoto system with equal edge weights, what is the largest possible link for a stable fixed point 
of the system? 


We consider the Kuramoto model as in (1.2 1 , where we assume that the graph is a ring with n vertices 
and that the edge weights are all equal, and by rescaling we can assume that they are all one, so that 
= 1 (here and below this sum on the indices is understood modulo n). Let 9* be a fixed point of ( |1.2[ ). 
From this we have 


0 = a;. + sin(0*+i - 9*) + sini9U - 0*), (4.1) 

and the Jacobian at this point is given by J, where 

=cos(0*-d4i), 

and the diagonal terms are chosen to make it a Laplacian. Then we see from the previous section that it is 
a stable fixed point if one of the following two conditions are met: 

1. For alH = I,..., n, 7 i,i+i := cos{9* - 9*_^_^) > 0; 

2 . one of the 7 qi+i is negative, and 


n—1 -| 

i =0 ''‘’*+1 




cos{9* - 9*^^) 


< 0 . 


(4.2) 


Note that cos(0) > 0 iff |d| < tt/2, and thus the sign of the cosine of the angle difference tells us whether 
a link is long or short. These conditions can be reinterpreted as saying: the fixed point is stable if all of the 
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links are short, or if there is one long link, then we need the sum of the reciprocals of the edge weights to be 
negative. The question then remains: is there a configuration 9* that is a fixed point of (1.2), has one long 
link, and satisfies (4.2)? 


Let us write Q = 0^+1 — 6^. Note that every term in the right-hand side of (4.1) will be a Q, except for 


the term 0 „_i — 0o which we can write as Ci- Then (4.2) becomes 

71 — 1 


E 


1 


cos(O) 


-f 


1 


cos 


(esc. 


< 0 . 


(4.3) 


Without loss of generality, let us assume that the longest link is given by the gap 0^-1 — 0o, so that if there 
is one long link, it is represented by the last term in (4.3). 

We can then pose this problem in an alternate form as follows: for any 77 G (0, 27r), let us find the 
maximum of 


n-2 


n—2 


/(C) = n 


2=0 


cos(Ci) 


subject to g{C) = ^ = 


V- 


i=0 


Then it is possible to realize a fixed point with largest gap 77 iff 

1 


mm/(C) 


cos 77 


< 0 . 


The nice thing we see from the constrained optimization problem is that since the constraint is symmetric 
up to any permutation of the variables, any internal optimizer must be as well. (We could also use Lagrange 
multipliers to obtain the same result.) Thus the maximum of / is acheived only when all of the Q are equal. 


and (4.3) becomes 


77—1 


cos(C) cos((77—1)C) 


< 0 . 


(4.4) 


Summarizing our conditions, if we assume that cos(C) is positive, then iff the conditions 

77—1 1 


cos ((77 — 1)C) > 0 or (cos ((77 — 1)C) < 0 and 


cos(C) cos((r7—1)C) 


< 0 


hold, then we have a stable fixed point for the solution 0^ = k(. Note that both of these conditions can be 
simplified to require that their product be positive, i.e. 




(4.5) 


To exhibit a “long link” solution, we need to show that there is C* with roots of hn{C*) = 0 ^nd \n(^* 
(mod 2)7r| G (7i'/2, 37r/2). We will show numerically that this is possible for every n, see Table 4.1 

Here we have found the first root of hn{C) and compute n(^n' this is the distance between 9q and 0„_i 
if we choose the solution 9k = kC,*. Note that in each case, the length of the link is greater than a quarter 
turn (note that we are presenting the length divided by 27r, so that a link is “long” iff it is greater than 1/4), 
but this distance decreases monotonically to 1/4 as 77 —)■ 00 . 

Note that in general, h„(C) has multiple roots, so that the first root may not always give the largest 
link, but we have found empirically that this is so. As such, the numbers in Table [4T] represent the longest 
link possible for a stable fixed point of the Kuramoto system on a ring of size n. 

We can also compute the values of ujk that give this solution. For all fc = 1,...,77 — 2, we have 
9k+i — 0fe = 0fc-i ~ = C) so that ojfc = 0. However, we also have 


a;„_i(C) := -sin(0o - 0„_i) - sin( 0„_2 - 0„_i) = sin(( 77 - 1)C) +sin(C), 


and Wo = —oJn-i- 
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n 

nC,* 12 'k 

3 

0.447 

4 

0.392 

5 

0.358 

10 

0.297 

20 

0.272 

30 

0.264 

40 

0.261 

50 

0.258 


Table 4.1 

Length of longest possible link for a stable solution on the ring with n vertices 



Fig. 4.1. We plot several of the functions referenced in the text for n = 9 vertices in the ring: 0.05 * hg{Cf) (blue), LosiO 
(orange), cos(8C) (purple). 


We plot several of the curves for the case of n = 9 in Figure 4.1 We have scaled / 19 (C) so that all the 


curves are of the same order of magnitude; note that only its sign matters. As shown above, where the blue 
curve is positive, the solution is stable. The purple curve is the cosine of the length of the (n — 1,0) link: 
since we see that this becomes negative before / 19 (C) does, there is a stable long link configuration. Finally, 
note that the zeros of / 19 (C) correspond to the critical points of a;8(C) (this can be checked), and thus the 
a;8(C) curve can be interpreted as a series of stable and unstable branches: whenever W8(C) is increasing in 
C, this corresponds to a stable solutions, and when it is decreasing, this corresponds to a 1-saddle. As we 
slide up and down in this graph, this leads to saddle-node bifurcations in the obvious manner. Moreover, 
note that the w = 0 slice is the one considered in 26 — it was shown there that for n = 9 there are two 


nontrivial sinks and two corresponding saddles, and this is recovered here. By decreasing or increasing wy, 
we can cause these saddles to collide with either of two sinks. 


5. Proof of Theorems. 

5.1. Proof of Theorem 12.61 

5.1.1. Definition of covering tree. Definition 5.1. Let G = (V, A, F) be a connected graph. We 
define the universal cover of G, denoted U{G). If G is a tree, then U(G) = G. If G is not a tree, then 
U{G) is a countably infinite tree defined as follows: 

• Choose a root vertex r G V(G). 
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• The vertices of U{G) are non-backtracking rooted paths in G, i.e. sequences of the form w = 
{r,V 2 ,V 3 ,... ,Vn) where (vt, Vi+i) G E, and v^-i Vi+i for all i. 

• We say that two vertices wi,W 2 G V{U{G)) are adjacent if one is an extension of the other, i.e. 

Wl = {r,V2,V3, . . .,Vn), W2 = {r,V2,V3, . . .,Vn,Vn+l), 

or vice versa. 

• Finally, we define the weight of the edge (wi,W 2 ) in the example above to be the same as the weight 
on the edge {vn-i,Vn). 

Definition 5.2 (Fundamental domain). For any graph G = {V,E,r), we define T, a fundamental 
domain of the cover of G, as follows. First construct the universal cover U(G). Choose T to be any finite 
connected subtree ofU{G) such that each edge of G appears once. Notice that by definition, G and T have 
the same edge sets, so we will write E for either set. Since there are \E\ edges in T and it is a tree, then 
\V[T) = \E\ + 1. We will also use the notation Nq = |y(G')| and Nt = |t^(T)| throughout. 

Remark 5.3. The universal cover has certain properties: 

1. If G is not a tree, then U{G) has countably many vertices. The degree of the node w = {r,V 2 , ■.. ,v„) 
in U{G) is the same as the degree of Vn in G. In particular, U{G) is locally finite. 

2. The map 

ip: U{G) G, 

{r,V 2 , ...,Vn)i-^Vn 

is a covering map, i.e. it is a map of graphs, it is surjective, and it is an isomorphism in the 
neighborhood of any vertex. 

3. The construction of T is independent (up to isomorphism) of the choice of the root vertex r. 

j. Finally, it should be noted that we can construct T directly without using the universal cover, see 
Example \5.5\ for an example. 

Proposition 5.4. Let G = (V,E,T) be a connected weighted graph, and define T as above. Then the 
graph map p: T ^ G given by 


p{{r,V2, . ■ ■,Vn)) = Vn 

is a surjective map. Note that |P(T)| — |1^(G)| = G, i.e. the difference in the number of vertices of a 
fundamental domain of the covering tree and the original graph itself is the cycle rank of the graph. 

Proof. Assume that wi,W 2 G V{T) and (wi,W 2 ) G ElT). Since wi,W 2 are adjacent in U{G), we have 

that 


Wl = (r,vi,.. .,Vk), W2 = {r,vi,.. .,Vk,Vk+i), 

and {vk,Vk+i) G ElJG). Then p{wi) = Vk and p(w 2 ) = vu+i are adjacent in G whenever wi,W 2 are adjacent 
in T. Therefore is a graph homomorphism. To see that it is surjective, note that every edge of G appears 
in T, and therefore its incident vertices must as well. Finally, note that 

\V{T)\-\ViG)\ = \E\ + l-\V{G)\=G. 


□ 

Example 5.5. Let us consider the diamond graph in Figure \5G\ If we choose the root vertex as 1, then 
2, 3,4 are neighbors and so we add 12,13,14 to T. Since 7;,, 7c have not yet appeared, we add them as in the 
tree. Of course, this is not the only covering tree we can obtain, even assuming that we choose 1 as the root; 
for example, if we remove vertex 143 and add a vertex 1234, this would also be a covering tree. 

5.1.2. The projections. Definition 5.6. Let G be a connected graph, and T, p: T ^ G as defined 
above. We define X to be the Nt x Nq matrix whose elements are 

Xa,v = l(<p(a) = v). 
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Fig. 5.1. The diamond graph G, and a choice of fundamental domain of its universal cover T. Note that G and T have 
the same edge sets with five edges, and T necessarily has six vertices. 


We will also refer to the columns by the vectors x'^^v G V{G). 

Lemma 5.7. Let G = {V^E^T) be a connected weighted graph, and T,X as in Definitions \5.2\\5.6\ Then 

Cg = X^CtX, (5.1) 

or 

{CG)v,^ = {x\CTxn- (5.2) 


Proof. Let us use formula (5.21. Writing this out, we have 


{x^Ctx^ =Y,xUCThxJ 
ij 


id 


If we can show that there exists a unique pair i,j such that (p{i) = v, ip{j) = w, and {CT)ij = then 

we are done. But notice that since ip{i) = v, and ip{j) = w, if {i,j) G E{T), then the weight on the edge 
{i,j) must be since is a graph homomorphism. Moreover, since each edge of G can appear at most 

once in T, this means that there is only one such pair. □ 

Definition 5.8. The map ip: T ^ G is surjective but in general not injective. For each v € V(G) with 
|(/7“^(u)| > 1, choose one representative of (p~^{v) as the primary representative of this set. Now, for 
each w G ip~^(v) with w prim(u), define the vector G by 


— 1, i = prim(ti), 
qf = I, i = w, 

0 , else. 


defined in Definition 5.6 Then the vectors q^ form a 


We define Q as the matrix whose columns are the q^ 

Lemma 5.9. Let S C R^^ be the span of the x'’ 
basis for S-^. In particular, dim(span(g“')) = rank((5) = c, the cycle rank of G. 

Proof. First we see that each q"^ is orthogonal to S. Consider x'" where v = <p{w). This vector is equal to 1 
both in the prim(i;) and w index, and therefore (g™, x"") = 0. Moreover, if we consider with v p{w), then 
this vector is equal to zero on these indices, and again {q^,x'") = 0. Since dim(span(5')) = rank(X) = |y(G)|, 
so that rank((3) = |i?| + 1 — dim S' = \E\ — |C(G)| + 1 = c. Finally, note that the g*" are clearly independent 
by construction, so by dimension counting they form a basis for S'*". 

□ 

Lemma 5.10. The Laplacian Ct'. > R^'^ is not invertible, but has a one-dimensional nullspace 

given by span(l). In this case, the canonical map 


Ct '. 


D At 


/{!} 


D At 


/{!} 
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is invertible. We then have that 


{Ct)-^Q = Zg. 


Proof. In general, for any graph H, the incidence matrix Bh is map from to In 

particular, since \E{T)\ = \E{G)\ and |I^(T’)| = \E{G)\ + 1, this means that 

is a map of co-rank one, and as such Bj. is a map with a one-dimensional nullspace, which is span(l). 
Therefore we can define B to be the canonical map so that 

is invertible: to define this, we first define B^ to be the induced map on the quotient and then B = (B^)^. 
In particular, note that B = B on 1-^. 

For each w such that ^ 0, consider the unique path in T from w to prim(i^(r(;)). Let us write this 
path as 

w = uo —>■ ui —>■ z ;2 Vn-i —>■ = prim(i^(w)). 

Define a vector r’" G as follows: 

{ 1 , e = {vi,v^+l),3i, 

- 1 , e = {vi+i,v^),3i, 

0, else. 

We first compute that 

B^r^ =q'". 

To see this, note that Bxr'^ = 0 on any vertex not in the path, and, moreover, will be zero on any vertex on 
the interior of the path due to cancellation. We need only consider the endpoints. 

First, assume that w ^ vi is the chosen orientation for the edge. Then {BTr'^)^ = ((J3 t)(w,vi)'''(w,vi) = 
1-1 = 1. If, on the other hand, vi w is the orientation, then we have {BTr'^)^ = (.{Bt)(w,vi)''’{w,vi) = 
-1 • -1 = 1 as well. Similarly, (i3T»’’")prim(vj(t«)) = “I- 
From this, it follows that if we define 

Y = B^^Q, (5.3) 


then the columi is of Y form a basis for the cycle space of G. 


By Lemma 2.5 Ct = —BTDBf, and it is easy to see that similarly, Ct = —BTDBf. From this we 


have 


giving the definition (2.1). 

□ 


Proof of Theorem 


2.6 


Q^{Ct)-^Q = -Q^B-^DB-^Q = -Y^ DY, 

Recall the definition of M := C(T). From Lemma 
Cg = X^CtX, 


5.7 


we have that 


and Lemma [A.3| gives that 


n+{CG) = n+{M\s), 
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since the columns of X form a basis for S. From Lemma [5.101 we have that 


Zg = Q^{Ct)-^Q, 


and Lemma |A.3| again gives that 


n+iZc) = n+{{M) 


since the columns of Q form a basis for S-^. Using Lemma A.l 


gives 


n+iCc) = n+iCr) - n+{ZG). 


Recall that T is a tree. Using Theorem 2.10 of [^, this implies that the number of positive eigenvalues 
of £t are given exactly by the number of negative edges in E(T). Recalling that the edge sets of G and T 
are the same, we are done. 

Finally, to prove the last sentence of the theorem, we again appeal to Theorem 2.10 of [^. One conse¬ 
quence of this theorem is that for £g to be stable, G+ must be connected. If the number of negative edges 
in G is the same as the number of cycles, this means that G+ is a tree, and thus is a spanning tree of G. 
For each negative edge in G, consider the two incident vertices and the unique path in G-|_ connecting them. 
This plus the addition of the negative edge make a cycle in G, and it is clear that each negative edge will be 
in a different cycle. Using this choice of cycles, we now have exactly one negative edge in each cycle. ■ 


5.2. Proof of Theorem |2.8[ We again start with several lemmas. 

Definition 5.11. Recall G,T,ip: T —>■ G. For each v G V(G), if = 1, let '4){v) = ip~^{v); 

else, 'tf{v) = prim(i;) as defined in Definition 5.8 Then order the vertices ofV{T) so that we list all of the 
vertices in the range of ip first, and then all of the other vertices in some order. It is easy to see that the 
range of ip as Nq elements. 

We then define the first Nq columns of the same way as we did in the definition of X, i.e. for v gV (G), 
define the vth column of U, u", by 


K)a 


1, g}{a)=v, 

0, else, 


and for columns i = Nq -|- 1,..., Nt, define m* to be the standard basis function with a 1 in the ith slot. By 
construction, U is a lower triangular matrix with Is on the diagonal, so det(17) = 1. 

Lemma 5.12. The Nq x Nq upper-left block of £tU is £g. 

Proof. The proof is more or less the same as Lemma [STT] Notice that once we have chosen this ordering, 
the first Nq columns of U is the matrix X as defined in Definition |5.6| with the corresponding ordering of 
vertices. □ 

Lemma 5.13. Let H be a symmetric matrix with a one-dimensional nullspace, and let x span this 
nullspace. Then for any y,z, we have 

det(H + i,»3)= 

{x,x) red 

Proof. Since H is self-adjoint, it has real eigenvalues and a orthonormal eigenbasis. Let us write Ai = 0 
with A 2 ,..., A„ not zero, and let wi,..., Wn be the corresponding eigenvectors. Let W be the matrix whose 
columns are the wp, since the eigenvectors form an orthonormal basis, W~^ = W^. 

Let Ml be the diagonal matrix whose diagonal entries are the Ai, and define M 2 = y) ® (W^z), i.e. 

(Majy = {wi,y) {wj,z) . 

Notice that H = WMiW^, since this is just its diagonalization, and 

WM 2 W^ = WiiW^u) 0 {W^w))W^ = WW^iu 0 w)WW^ =u(^w, 
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where we have used the identity (Ax) 0 (By) = A{x ® y)B"^. Thus we have 

H + u®w = W{Mi + M 2 )W^, 


and since W is unitary this gives 


det{H + u(S>w) = det(Mi + M 2 ). 

To compute the last term, notice that we can use the multilinearity of the determinant by rows and break 
up this into 2" terms. However, note two things: the matrix M 2 is rank one and thus has only one linearly 
independent row, so if we choose more than one row from M 2 , the determinant is zero. Also, the first row 
of Ml is zero, so to obtain a nonzero determinant, we must choose the first row of M 2 . Therefore the only 
term which gives a nonzero determinant is the one where we choose the first row from M 2 and all of the 
other rows from Mi. Expanding by minors in the first column, the determinant of this matrix is 


{wi,y) {wi,z) A 2 X • • • X A„. 

Since we chose Ai = 0, this means that Wi is the unit vector that spans ker(M), and also notice that detred H 
is by definition A 2 x • • • x A„. Thus we have 

det(M + y®z) = {wi,y) {wi,z) det H 

red 

If a; is a general vector in ker(i/), then x is a scalar multiple of wi and we have 

det(H + i,»3)= 

{x,x) red 


□ 

Proof of Theorem 12.81 Recall the definition of U in Definition 15.111 Define 
A = U^{Ct + 1<S) 1)U = U^CtU + U^l (g) U^l. 

If we can establish the following three identities: 

det (A) = 7 e, 

eeE(G) 

det(A|s) = det Cg, 

I\Q red 

det((A"i)|s^) = det(ZG), 


then, using Lemma |A.2[ we obtain 




n >= 

eeE(G) 


—— QStred 

I\g 


det (Zg) 

and this proves the theorem. So it remains to show these three identities. 


Identity (5.6) follows directly from Lemma 5.10 To establish (5.4), we note that 
det(A) = det{U ^{£t + 1 0 d)U) = det(£T + 10 1), 


since det(C/) = 1. By Lemma 5.13[ we have 

det(£T 


1 (g) 1) = ^ ^ ^ det Tj' = A^t’ det £tj 

( 1 , 1 ) red red 
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and since there is only one spanning tree of T, by the Matrix Tree Theorem, 


detCT = NT TT 7e- 

red 

e£E{G) 


To establish (5.5), note that if we use the basis for made from the columns of f/, then A\s is ju st the 
upper left Nq x Nq block of A, which is itself the upper left block of CtU +By Lemma 5.12 


the upper left block of the first term is just Cqj and the upper left block of the second term can be written 
as 




where ([/) is just the first Nq columns of U. In particular, the vth row of t/^1 is a zero-one vector with as 
many ones as vertices in and therefore u := is a column vector of height Nq whose nth entry 

is |<^"^(u)|. 

Putting this together, and again using Lemma [5. 13| gives 


det(A| 5 ) = det(£G +u®u) 


{u, 1) {u, 1) 

( 1 , 1 ) 


det Cg 

red 


(E.gV(G)£ 

Ng 


det/lc = 

red 


iV| 


det Cg ■ 

red 


Appendix A. Classical lemmae. 

This follows from the following duality formula for the number of negative eigenvalues (or more generally 
the inertia) of a matrix. This result seems to have been frequently rediscovered but the first instance we 
have found in the literature is due to Haynsworth 31 . This result has been used frequently in the nonlinear 


waves literature, in the context of the stability of traveling wave solutions. 

Lemma A.l (Haynsworth). Suppose that M is a non-singular NxN Hermitian matrix. Let S a subspace 
of subspace ofM.^, S'^ the orthogonal subspace, Ps the orthogonal projection onto S, and M\g = PgMPg 
the restriction of M to S. If M\g is non-singular, then 

n+{M) = n+{M\s) + n+{{M-^)\s±) 


Note that, in her original paper, Haynsworth states this theorem slightly differently in terms of the Schur 
complement of the matrix, but it is obviously equivalent. 

Lemma A.2. Suppose that M is a non-singular N x N Hermitian matrix. Let S a subspace of sub space 
ofM.^, S'^ the orthogonal subspace, Pg the orthogonal projection onto S, and M\g = PgMPg the restriction 
of M to S. If M\g is non-singular, then 


det(M) 


det(M|s) 

det(M-i|si)' 


This is more or less the determinantal version of Lemma and could be proved in a similar fashion. 
It can also be found explicitly in . 

Lemma A. 3 (Sylvester’s Law of Inertia). Suppose M is Hermitian and U is non-singular. Then 

n_(M) = n_{U^MU). (A.l) 

In particular, if A is an n x n matrix, S C M” is a k-dimensional subspace, Pg is the orthogonal projection 
onto S, xi,... ,Xk is a basis for S, and X = {xi, ..., Xk) is the k x k matrix whose ith column is Xi, then 

n+iPgAPg)=n+iX^AX). 

Proof. The first statement is just the standard formulation of Sylvester’s Theorem, so we will prove the 
second part. Notice that if we choose y^,... ,y^ to be an orthonormal basis for S, then PgAPg = AY. 
Choosing Z = Y~^X, we have 

Z^{Y^AY)Z = X^AX, 

and by the first statement this implies that these two matrices have the same signature. □ 
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